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Abstract 

We analyze the use of layered superconductors as anisotropic metamaterials. Layered super- 
conductors can have a negative refraction index in a wide frequency range for arbitrary incident 
angles. Indeed, low-T c (s-wave) superconductors allow to produce artificial heterostructures with 
low losses for T <C T c . However, the real part of their in-plane effective permittivity is very large. 
Moreover, even at low temperatures, layered high-T c superconductors have a large in-plane normal 
conductivity, producing large losses (due to d-wave symmetry). Therefore, it is difficult to enhance 
the evanescent modes in either low-T c or high-T c superconductors. 

PACS numbers: 74.25.Nf, 42.25.Bs 
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Metamaterials are attracting considerable attention because of their unusual interaction 
with electromagnetic waves (see e.g., [l]). In particular, metamaterials supporting nega- 
tive refractive index have the potential for subwavelength resolution 2] and aberration-free 
imaging. The first proposed negative index metamaterials used subwavelength electric and 
magnetic structures to achieve simultaneously negative permittivity e and permeability fi 
(see, e.g., ^]). However, these "double negative" structures require intricate design and 
demanding fabrication techniques, are not "very subwavelength", and suffer from spatial 
dispersion effects. Moreover, the implicit overlapping electric and magnetic resonances (see, 
e.g., Q]) often leads to resonant losses that, together with material losses, lead to significant 
degradation in metamaterial functionality. This manifestation of loss can be quantified by 
examining the figure of merit (FOM) in such materials, which is defined as \n'\/n" where n' 
and n" are the real and imaginary parts of the refractive index n, respectively. The FOMs 
of negative index materials in the visible and near-IR have experimentally ranged from 0.1 
up to 3.5 flfl. 

Another promising route to creating negative index metamaterials is to use strongly 
anisotropic materials, in particular, uniaxial anisotropic materials with different signs of the 
permittivity tensor components along, en, and transverse, e±, to the surface (see, e.g., 
These materials have been theoretically pi and experimentally [9] demonstrated to support 
sub- wavelength imaging, and they have also been proposed as a model system for scattering- 



free plasmonic optics 



lOj and subwavelength-scale waveguiding [11]. These materials are 



particularly attractive because: they are relatively straightforward to fabricate, compared 
to double negative metamaterials; they do not require negative permeability; and do not 
suffer from magnetic resonance losses. The FOMs for such materials have been calculated 
to be significantly greater than those measured in double negative materials 

Experimental schemes for creating strongly anisotropic uniaxial materials have typically 
involved the fabrication of subwavelength stacks of materials whose layers comprise alternat- 
ing signs of permittivity. For example, alternating stacks of Ag and AI2O3 9] and of doped 
and undoped semiconductors [12] have been demonstrated to support strong anisotropy 
in the visible and infrared frequency ranges respectively. However, spatial dispersion can 
strongly modify the optical response of the system relative to the ideal effective medium limit 
response (l3j; strong local field variations exist due to the structure and length scales of plas- 
monic modes supported by negative-permittivity films, even in the limit of I <C clq, where I 



2 



is the length scale of the thin films in the material and ao is the free-space electromagnetic 
wavelength. This imposes limitations to subwavelength imaging and waveguiding in such 
materials. Spatial dispersion may be reduced by making composite structures with thinner 
layers. However, there exist practical material deposition limitations to thin-film stacks in- 
volving film roughness and continuity. In addition, damping due to electron scattering at 
the thin film interface becomes significant starting at length scales of aoVp/c ~ ao/100 where 
v f is the Fermi velocity in the material [ijj] , limiting the minimum film thickness. It is clear 
that composite structures are limited in practice as "ideal" strongly anisotropic materials. 
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We analyze here the idea of using superconductors as metamaterials (see, e.g., [15|, 



17]). In particular, we consider layered cuprate superconductors [17| and artificial 
superconducting-insulator systems [3] as candidates for strongly anisotropic metamate- 
rials. Unlike the composite structures discussed earlier, layered superconductors are not 
limited in performance by the spatial dispersion effects discussed in 13j. We will analyze 
these materials in the specific context of subwavelength resolution, which can be achieved 
by the amplification of evanescent waves [2| . This amplification is high when n is close to 
unity and its imaginary part is small j^, 19]. For the incident p-polarized waves considered 
here, subwavelength resolution requires Im(e) <C exp(— 2kj_L), where k± is the wavevector 
component across the surface, and L is the plane lens-thickness 19J. For evanescent modes 
with kj_ = 2u/c = 2ko and L/a^ = 0.1, we have Im(e) -C 0.081. 

We show that in the case of natural high-T c cuprates the losses are high at any reasonable 
frequency. In the case of artificial layered structures prepared from low-T c superconductors, 
the losses can be reduced significantly at low temperatures, T -C T c , where T c is the critical 
temperature. The frequency range for such a metamaterial is hu < 2A, where A is the 
superconducting gap, which corresponds to a maximum frequency in the THz range for 
low-T c superconductors. We prove that the in-plane permittivity for low-T c multi-layers 
is large, preventing the effective enhancement of evanescent waves. This is problematic 
because subwavelength resolution 2j requires the amplification of evanescent waves. Note 



that Refs. 
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only focus on the zero-frequency DC case. 
Effective permittivity. — We study a medium consisting of a periodic stack of supercon- 
ducting layers of thickness s and insulating layers of thickness d with Josephson coupling be- 
tween successive superconducting planes. The number of layers is large, L/(s + d) = N ^> 1 
and s is smaller than: the in-plane magnetic field penetration depth An, transverse skin 
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depth 5±(u), and wavelength a(u) ~ 2-Kc/uyJ\e{u)\. We calculate the effective permittivity, 
e = (£ii,£_l), of the layered system in the case of p-wave refraction. 

Layered superconductors with Josephson couplings can be described by the Lawrence- 
Doniach model, where the averaged current components can be expressed as {^(J 

J ^ = JcSm ^ + 2^(7Trf)' J|! = 8^AJ +CT||E|1 ' (1) 
where y? n is the gauge- invariant phase difference between the (n + l)th and nth supercon- 
ducting layers, p n is the in-plane superconducting momentum, J c = c$o/(8vr 2 (iA 2 L ) is the 
transverse supercurrent density, $o is the magnetic flux quantum, and Aj_ is the transverse 
magnetic field penetration depths. Also a± and an are the averaged transverse and in-plane 
quasiparticle conductivities. The transverse E± and in-plane E\\ components of the electric 
field are related to the gauge-invariant phase difference and superconducting momentum 



by |2fl|,|2l| 

where V^/(n) = f(n + l) + f(n — 1) — 2f(n), a = eR 2 D /(sd) is the capacitive coupling 
between layers, and Rd is the Debye length. We linearize the first of Eqs. ([1]) and consider 
a linear electromagnetic wave E^±(x,n,t) = ^J'j^zE^^k, q, to) exp(— iut + ikx + iqn), 
where q = nl/(N + 1), I = 0,±1,±2, and the x-axis is in the plane of the layers. Using 
Eqs. ([I]) and (j2J), we obtain: 



^=(l + aq* 



eu, 



a x - 



< 2 p (s+d) 



2, ,2 



J ll £ T u p 



4iridu 

where u> p = cj (\j_y/e) is the Josephson plasma frequency, e is the interlayer permit- 
tivity, 7 = Aj_/A||, and q 2 = 2(1 — cosg). Averaged over the sample volume, the 
Maxwell equation has the form cV x H = Arc J + dD/dt, where D\\ = b2E\\ and 
D±_ = e° ± E±. In the effective medium approximation, the components of the permittiv- 
ity tensor can be expressed as [22| ejj = (de + s)/(s + d), and £° = e(s + d)/(s£ + d), 
where we assume that £ superconductor = 1- Fourier transforming the above Maxwell equa- 
tion, we derive c [V x H] ± (k, q, uj) = —ej_Ej_ and c [V x Hi, (k, q, u) = —e\\E\\, where 
£|| = ejj — (&ir/iu)(J\\/E\\) and £± = £± — (4:ir/iuj)(J±/ E±) . Therefore, we finally obtain 



e± = e± : e(l + aq . 

too V / u z d 

Airan „<Jl 



Thus, £ii < and e± > in the frequency range 



~\se + d 00 e(s + d) , N 

1 + aq 2 ) — — < — <7\ M , • 5 
/ a uj p V as + s 

If the incident angle is close to normal and anisotropy is large, 7 > 1, we can find an estimate 

FOM « 2 Re(e||)/Im(e||)| ~ £7 2 a; 3 /27rcr||u; 2 . Electromagnetic waves propagate in the layered 

superconductors if uj > u p . Thus, the results obtained are valid if u p < uj < u c = 2A/h. 

Below we analyze separately the different cases of a typical high-T c layered superconductor, 

E^S^CaC^Os+s (Bi2212), and also an artificial low-T c layered structure made from Nb. 

Layered high-T c superconductors. — In the case of Bi2212, it is known that s < (I = 1- 

2 nm, e = 12, a w 0.1, and at low temperatures (T < T\ = 90 K) u p » 10 12 s _1 , 7 = 500, 



« 4 • 10 4 ft^cm -1 , and ct_l « 2 • 10~ 3 fi^cm" 1 [20|, [23J. In this case, Eqs. flU can 
be rewritten as e± we(l- u; 2 /u; 2 ) +47T20"j_/u;, ey 7 2 u; 2 /u; 2 ) +47rzo"|j/u;. The calculated 

frequency dependence of the permittivity for Bi2212 is shown in Fig. 1. The superconducting 
gap for Bi2212 is estimated as A m 2-3/cb^c, with u c ~ 5 x 10 13 s _1 7^. Thus, for 
any incident angle, Bi2212 has negative n in the frequency range from about 0.15 THz to 
7.5 THz, or in the wavelength domain 40 /im< a < 2 mm. However, the use of Bi2212 as 
metamaterial has a disadvantage since the in-plane quasiparticle conductivity cni is large, 
even at helium temperatures. As it is seen from the inset in Fig. 1(b), o~\\ 7^ when T — ► 0, 
which is typical for superconductors having a d-type symmetry of the order parameter. In 
addition, the usual dimensions of high-quality Bi2212 single crystals are less than 1 mm in 
the in-plane direction and about 30-100 /im in the transverse direction. Thus, it might be 
difficult to use Bi2212 single crystals as metamaterials, or elements of a superlens. 

Low-T c artificial layered structures. — The thickness of the insulator in Josephson junc- 
tions is about a few nm. To attain a low-loss regime and reach the bulk critical temper- 
ature, the thickness of the superconducting layers should be larger or about the super- 
conductor coherence length £. For clean superconductors, £ is about tens of nm. Thus, 
for low-T c artificial-layered structures, it is reasonable to analyze the case d s. In this 
limit, An = A-v/(s + d)/s ~ A, where A is the bulk magnetic field penetration depth and 
a = eR 2 D /(sd) < 1 in any realistic case. It is easy to choose an insulator with very 
low conductivity o~i to satisfy the condition crj o~ s d/s at any reasonable temperature, 
where a s is the quasiparticle conductivity of the superconductor. In this case we have 
e° ± = 1, e?i = 1, cr_i_ = Oisjd and an = cr s . Equations (j3J for the effective permittivity can 
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FIG. 1: Dependence of the real and imaginary parts of the permittivity e in Bi2212 on frequency 
u (or wavelength eto), calculated from Eqs. (j3|): (a) real part of the in-plane permittivity £n(w). 
Inset: ratio of the real parts of the in-plane and transverse permittivities; (b) imaginary part of the 
in-plane permittivity. Inset: temperature dependence of the in-plane quasiparticle conductivity 
o"ii(T) in Bi2212; solid tr iang les are low frequency data from Ref. |23| . open squares correspond to 



14.4 GHz data from Ref. |24; (c) imaginary part of the transverse permittivity £±(lo). 
now be rewritten as 

Therefore, the refraction index n is negative if 

\Jes/d < cu/ujp < 7. (7) 

For artificial structures, 7 can be easily made of the order of, or even much larger than, in 
natural layered superconductors. In contrast to c?-wave high-T c superconductors, for bulk 
s-wave superconductors, the quasiparticle conductivity a s tends to zero for decreasing T. 




FIG. 2: Calculated, from Eq. (jHJ), temperature dependence of the imaginary part of £\\(t) = 
£n(T/T c ), for a Nb-based layered structure, with to = 0.9u> c , e = 10, s/d = 5, and 7 = 500; here: 



lo p /lo c = 0.1, Re (e±) = 0.393, and Re (e» 



-3 • 10 6 . The inset shows the dependence 



25] of 



(J s ja n on i e T/T c ; points: experimental data for Nb at about 60 GHz; solid line: Mattis-Bardeen 



theory in the weak-coupling BCS limit; dashed line: strong-coupling Eliashberg prediction 



25|]. 



Thus, in principle, the imaginary part of £11 could be made as small as necessary by cooling 
the system. 



Consider now Nb superconducting layers. For estimates we can take 25j: T c = 9.3 K, 
A(T = 0) = 44 nm, £ = 38 nm, electron mean free path l e = 20 nm, and normal state 
conductivity a n = 0.85 x 10 6 f2 -1 cm -1 . Thus, a reasonable thickness for the superconducting 
layers can be chosen as s = 30-40 nm -C a(u c ) ^ 100-200 nm. Superconducting properties 



n particular, its 
2611 (see inset in 



of Nb are well described in the BCS weak-coupling approximation [251 ] . 
conductivity a s (u,T) can be calculated using the Mattis-Bardeen theory 
Fig. 2). At low temperatures, T -C T c , in the weak-coupling BCS limit, we have A 
1.76 ksT c . When u < u c and T T c , we can rewrite the Mattis-Bardeen formula for 
conductivity 



25 



261 ] in the form 



Vs/on = [1 — ex P (— 3.52cj/ (uj c t))] /uj X 



00 (u 2 



(w 2 - 1) [(u + 2uj/uj c ) 2 - 1] 



du, (8) 



where t = T/T c . The results of our calculations are shown in Fig. [21 These calculations 
demonstrate that the losses in artificial structures made from low-T c superconductors can 
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be extremely low. The maximum frequency u c = 3.52 ksTc/h for Nb corresponds to approx- 
imately 0.7 THz. From the results presented in Fig. 2, we can estimate that at uj ~ uj c 
the imaginary part of en is lower than 10 -3 if T < 1 K. At higher frequencies, to > u c , 
the conductivity of the superconductor is about the conductivity of the normal metal and 
it cannot be easily used as a metamaterial with low losses. Note also that by an ap- 
propriate choice of insulator, s, and d, we can vary the parameters 7 and u p in a wide 
range. If we assume that s ~ 10, then to fulfill conditions (jTJ) for uj p < uj c we should 
prepare highly-anisotropic heterostructures with 7 > 10 3 . If the anisotropy is large, we 
can find from Eq. ([6]) that Re (e\\) ~ —c 2 /\ 2 uj 2 . The absolute value of Re (e\\) is very 
large, |Re (eii)| > c 2 /\ 2 oj 2 ~3x 10 6 . These estimates suggest that low-T c superconducting 
multi-layers might not work as practical metamaterials. 

The metamaterial properties of layered superconductors, either natural or artificial, can 
be tuned varying the temperature or an in-plane magnetic field, which strongly affects the 
transverse critical current density and, consequently, the plasma frequency. But applying 
a magnetic field increases dissipation, which is undesirable. Note also that the estimates 
made above show that the value of FOM may be very large for the systems considered here, 
however, this does not mean necessarily that these media can be easily used as practical 
metamaterials. 

Cuprates in the normal state. — There is experimental evidence that cuprate supercon- 
ductors have strongly anisotropic optical characteristics in the normal state 271, 128]. For 
example, it was observed that L^-zSr^CuC^ supports negative permittivity along the CuO 
planes at frequencies up to the mid- and near-IR range 27j . Moreover, these optical prop- 
erties could be finely tuned by varying the stoichiometry. Such natural materials are thus 
candidates for practical anisotropic metamaterials. The use of cuprates in the normal state 
have evident advantages, such as operating above u c and to work at room temperature. 
However, the normal conductivity of cuprates is of the same order as their quasi-particle 
conductivity in the superconducting state (see, e.g., the inset in Fig. lb and Ref. |23j). The 
metamaterial properties of cuprates in the normal state require a separate analysis and will 
be performed elsewhere. 

Conclusions. — Here we analyze the properties of anisotropic metamaterials made from 
layered superconductors. We show that these materials can have a negative refraction index 
in a wide frequency range for arbitrary incident angles. However, superconducting metama- 
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terials made from natural layered high-T c cuprates have a large in-plane normal conductivity, 
even at very low temperatures, due to rf-wave symmetry of their superconducting order pa- 
rameter. Therefore, these are very lossy. Nevertheless, low-T c s-wave superconductors allow 
to produce metamaterials with low losses at low temperatures, T <C T c . But the real part of 
their in-plane permittivity is very large, reducing the enhancement of the evanescent modes 
and potentially limiting the use of superconducting structures as practical metamaterials. 

We gratefully acknowledge partial support from the NSA, LPS, ARO, NSF grant No. 
EIA-0130383, and JSPS-RFBR 09-02-92114, FC gratefully acknowledges useful discussions 
with E. Narimanov. 
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